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A decentralized LQG-based control strategy is designed to reduce low-frequency sound transmission through 
periodically stiffened panels. While modem control strategies have been used to reduce sound radiation from 
relatively simple structural acoustic systems, significant implementation issues have to be addressed before these 
control strategies can be extended to large systems such as the fuselage of an aircraft. For instance, centralized 
approaches typically require a high level of connectivity and are computationally intensive, while decentralized 
strategies face stability problems caused by the unmodeled interaction between neighboring control units. Since 
accurate uncertainty bounds are not known a priori , it is difficult to ensure the decentralized control system will 
be robust without making the controller overly conservative. Therefore an iterative approach is suggested, which 
utilizes frequency- shaped loop recovery. The approach accounts for modeling error introduced by neighboring 
control loops, requires no communication between subsystems, and is relatively simple. The control strategy is 
validated using real-time control experiments performed on a built-up aluminum test structure representative of 
the fuselage of an aircraft. Experiments demonstrate that the iterative approach is capable of achieving 12 dB peak 
reductions and a 3.6 dB integrated reduction in radiated sound power from the stiffened panel. 


1 Introduction 


2 Decentralized Control 


Interior noise levels in commercial and general aviation air- 
craft are an occupational hazard for flight crew and an incon- 
venience for passengers. Since passive treatments, such as 
poroelastic foam, are not effective at low frequencies due to 
size and weight constraints, this paper focuses on an active 
control strategy. In particular, primary structural control is 
considered. As the name implies, this approach reduces inte- 
rior noise by applying control inputs directly to the primary 
structure (e.g., piezoelectric patches bonded to the aircraft 
sidewall). 

Past research has shown that active structural control can be 
used effectively on laboratory scale systems. However, sig- 
nificant implementation issues have to be addressed before 
active control can be implemented on a full-scale aircraft 
[ 1 ]. While centralized strategies have been used to simul- 
taneously reduce the sound radiation from multiple aircraft- 
style bays [ 2 ], these strategies require a high level of con- 
nectivity, are computationally intensive, and can be sensi- 
tive to transducer failures. On the other hand decentralized 
strategies are simple, computationally efficient, and scalable. 
Unfortunately the performance of decentralized strategies is 
often limited due to the interaction between neighboring con- 
trollers [ 2 ] . The goal of this work is to show that a scalable 
decentralized control system can be effective when imple- 
mented in a simple iterative loop. 

This paper begins with a brief description of decentralized 
control along with a discussion of its limitations. In particu- 
lar, we demonstrate that the approach is only effective if the 
local control units are robust to the modeling error introduced 
by neighboring controllers. While it might be tempting to 
use inherently robust strategies like direct velocity feedback, 
the conservative nature of these approaches provides stabil- 
ity guarantees at the expense of performance. Therefore this 
paper focuses on linear quadratic Gaussian (LQG) control. 
Since LQG designs can have arbitrarily poor stability mar- 
gins, frequency shaping and loop transfer recovery (LTR) 
are often necessary to improve the robust stability and per- 
formance of the system. However, both of these techniques 
assume uncertainty bounds are known a priori , which is not 
always true. Therefore an iterative scheme is presented to ac- 
count for the lack of a priori information. In addition to dis- 
cussing the control methodology, this paper also describes 
the experiment that was used to evaluate the proposed ap- 
proach. 


Decentralized control implies that each control unit is de- 
signed and implemented independently, as demonstrated by 
Fig. 1 (a). In this example the plant is defined as 


G = 


Gn G 12 

G21 G22 


( 1 ) 


where Gij denotes the transfer function from the it h actua- 
tor to the jth sensor. In the context of this work, the diago- 
nal terms of G model the response between actuator/sensor 
pairs on the same bay, while the off-diagonal terms capture 
the coupling between actuators and sensors on different bays. 
Decentralized control is particularly effective if the plant is 
diagonally dominant (ie G11G22 » G12G21) since the 
cross-coupling between the ith input and the jth output is 
neglected during the design process. Figure 1 (a) can also be 
rearranged as shown in Fig. 1 (b). Notice that while the con- 
troller C2 is designed for the nominal subsystem model G22 , 
the cross-coupling terms (G 12 and G21) and the controller 
Ci introduce an additional path from U2 to 7/2- Therefore, 
the combined model from U2 to 7/2 is 


G2 — G22 + ( 2 ) 


where A^i = G21C1G12/ (1 — GiGn) is the additive er- 
ror introduced by C%. If the local control system designed 
for G22 is not robust to the modeling error introduced by Gi, 
then the coupled system will be unstable. Therefore, decen- 
tralized control requires robust local controllers. 

Robust stability can be expressed in terms of the gain margin, 
which is a measure of the system’s tolerance to multiplicative 
error. Therefore, it is beneficial to rewrite Eq ( 2 ) as 

G2 = G22 (1 + Ami) ( 3 ) 


where (1 + Ami) is the multiplicative error, and the multi- 
plicative uncertainty Ami is defined as Ami = A^i/G 22- 
In this example, 


Ami 


(G 2 iG 12 \ f GxGn \ 
VG11G22; \1 — GiGn / 


( 4 ) 


Notice that the first term in parentheses provides a measure 
of the diagonal dominance of the plant. This term will be 
large if the plant is not diagonally dominant (ie the product 
of the cross-coupling terms is large with respect to the prod- 
uct of the diagonal terms). The second term in Eq ( 4 ) is the 
complementary sensitivity function, T \ . To maintain robust 




Figure 1 : (a) Standard block diagram for a two-channel 
decentralized control system; (b) Rearranged diagram 
highlighting the additive error. 


stability, the control system designed for subsystem 2 must 
be robust in the frequency bands where the plant is not diag- 
onally dominant and at frequencies where T\ is large. Since 
decentralized control units are designed completely indepen- 
dently, there is no way of knowing the cross terms or the 
complementary sensitivity function in advance. However, we 
do have access to the local control models, G\\ and £ 22 - No- 
tice in Eq (4) that the multiplicative uncertainty is inversely 
proportional to the product of the diagonal terms. Therefore, 
the uncertainty tends to be large at the zeros of Gn and 622 - 
As a result, it is often advantageous to penalize control effort 
at the zeros of the local control model. The following section 
explains how a frequency- shaped LQG design can be used to 
achieve this goal. 


3 Linear quadratic Gaussian (LQG) 
control 

LQG control is an optimal control strategy that uses mini- 
mum variance state estimates with an optimal state regulator. 
While the state regulator has guaranteed stability margins, 
the LQG controller can have arbitrarily poor stability margins 
due to errors in the state estimates [3]. Therefore, frequency 
shaping and loop transfer recovery are both commonly used 
to improve the performance and stability of the design. 

Frequency- shaped designs can be generated by shaping the 
state regulator, the stochastic estimator, or both [4] . For in- 
stance, one way to shape the regulator is to augment the plant 
model with filter dynamics that produce additional frequency- 
weighted outputs. The frequency- shaped outputs are then in- 
cluded in the cost function and the optimal feedback gain 
matrix is calculated to minimize this new function. There- 
fore, we can penalize the control effort at the zeros of the 
local control model by augmenting the plant with low-order 
band-pass filters that generate an additional output Uf. The 
cost function can then be expressed in terms of this new out- 
put as 


J = 



(■ q\y 2 + q 2 U 2 f + ru 2 ) dt 


( 5 ) 


where q\ is the original output weighting term, q<± weights the 


new output, and r weights the frequency independent con- 
trol effort u. In this case, the effective effort penalty at each 
frequency will be the sum of the frequency dependent and 
independent terms [5]. 

While it is often useful to penalize control effort using fre- 
quency shaped designs, modeling error is not restricted to 
the zeros of the local control models. Therefore it is not al- 
ways possible to penalize control effort and improve robust- 
ness without significant performance penalties. Loop trans- 
fer recovery provides an alternative that can often be used 
to improve robustness without making the controller overly 
conservative. 


4 Loop transfer recovery (LTR) 

As discussed earlier, the guaranteed stability margins asso- 
ciated with full- state feedback are lost when the estimator 
is introduced. In response to this problem, loop transfer re- 
covery (LTR) was developed by Doyle [ 6 ] to asymptotically 
’’recover” the properties of the full-state feedback design. Es- 
sentially, LTR modifies the estimator such that the LQG sys- 
tem inherits the robustness of the state regulator. As the name 
implies, the robustness of the regulator is recovered by mak- 
ing the LQG loop transfer function approach that of the state 
regulator [7]. Recovery is achieved by injecting fictitious 
noise at the plant input. If the plant is stabilizable, completely 
observable, time-invariant and also minimum phase, then the 
LQG design will asymptotically recover the characteristics 
of the state regulator as the amplitude of the fictitious input 
noise becomes arbitrarily large [4] . 

Unfortunately, full loop transfer recovery is only guaranteed 
if the system is minimum phase. This is because full recovery 
relies on pole-zero cancelations. Note that even if the plant 
G = C(sl — A)~ 1 B is not minimum phase, the continuous- 
time LQR loop transfer function, L lqr = — K (si — A)~ 1 B 
is minimum phase. This explains the attractive gain and 
phase margins associated with the full state design. However, 
if G is not minimum phase, then full recovery would im- 
ply that the compensator uses right-half-plane poles to cancel 
the right-half-plane zeros. Since no stabilizing compensator 
can use unstable poles to cancel right-half-plane zeros, full 
recovery is only possible if the system is minimum phase. 
Unfortunately, noncollocated transducers often yield right- 
half plane zeros. Therefore, full recovery is rarely possi- 
ble. Luckily, partial loop recovery is often adequate to ac- 
count for plant uncertainty. Partial recovery is achieved by 
increasing the amplitude of the fictitious input noise until 
maximum robustness is achieved. Another option is to use 
frequency-shaped recovery. This enables different perfor- 
mance/robustness trade-offs in different frequency bands [ 8 ]. 

Although loop recovery can be effective, it is limited to re- 
covering the characteristics of the state regulator. Therefore, 
the approach is only suitable if the underlying LQR design 
is robust. As a result, a two-step approach is advocated. 
First, the effort penalty is shaped to ensure the LQR design 
is sufficiently robust. This is accomplished by limiting con- 
trol authority at the zeros of the local control model. As 
discussed earlier, this typically improves the robustness of 
the design without significant reductions in achievable per- 
formance. Next, frequency- shaped loop recovery is used in 
frequency bands where the interaction between controllers 











Figure 2: A diagram of the iterative control strategy. 


limits achievable performance. Unfortunately, it is often dif- 
ficult to accurately predict the frequency bands where loop 
recovery is required a priori. Therefore the following sec- 
tion describes the concept of iterative loop recovery. 

5 Iterative loop recovery 

The iterative approach described in this paper incorporates 
loop recovery in a simple iterative loop. The strategy is de- 
picted in Fig. 2 using Gi(z) to represent the physical system 
and di (n) as the disturbance. The subscript i denotes that the 
control strategy is implemented on the i th subsystem. The 
controller is designed using an internal plant model Gi(z ), 
which is used to generate an estimate of the disturbance di (n) 
at each time step. Although this resembles an internal model 
control structure [9], the input to the online controller Ci(z) 
is the observed error signal yi(n) instead of the disturbance 
estimate. While the diagram only depicts one control unit, 
the control strategy is implemented simultaneously on each 
subsystem. 

The control strategy can be divided into two processes: a 
real-time process that occurs on a sample-by-sample basis, 
and an update procedure that occurs much less frequently. 
The steps involved in the update procedure are connected by 
the dashed lines in Fig. 2. Essentially, an LQG controller 
is updated periodically using a new effort weighting term 
and disturbance model. The disturbance model is calculated 
from the discrete Fourier transform of the disturbance esti- 
mate, Di(k). Once the disturbance model and effort weight- 
ing term are updated, the LQG controller is redesigned. If 
the design satisfies control effort and performance consid- 
erations, then the online controller Ci{z) is updated. The 
entire process is repeated until the control effort approaches 
predefined limits, or until subsequent designs fail to improve 
control performance. 

While it might not be obvious at first, this strategy relies on 
frequency shaped loop recovery. Loop recovery is exploited 
by increasing the amplitude of the disturbance model. Al- 


though changes in the disturbance model introduce fictitious 
noise at the plant output, output noise can also be represented 
as appropriately shaped input noise if the system only has 
one input and one output. Therefore, increasing the ampli- 
tude of the disturbance model is used to improve the robust 
stability and performance of the control system by partially 
recovering the characteristics of the LQR design. 

Recall that the proposed control strategy uses an internal 
model to generate the disturbance estimate di(ri). If the 
plant model is perfect, then the estimate will accurately track 
changes in the true disturbance. However in frequency bands 
where Gi(z) ^ Gi{z ), the amplitude of the disturbance es- 
timate could be larger or smaller than the actual disturbance. 
While increasing the amplitude of the disturbance model typ- 
ically improves robustness by exploiting loop recovery, de- 
creasing the amplitude of the model can make the design 
more sensitive to modeling error. Therefore the disturbance 
model is updated by taking the maximum of the current dis- 
turbance estimate and the previous disturbance model on a 
frequency-by-frequency basis. In other words, the model is 
only changed if the magnitude of the disturbance estimate 
exceeds the magnitude of the previous disturbance model. 

In summary, the first step in the design process is to iden- 
tify an accurate subsystem model. The effort weighting term 
is then shaped to penalize control effort at the zeros of the 
open-loop system. Shaping the effort term limits excessive 
compensator gain and improves the stability margins of the 
LQR design. Next, the measured response is used to generate 
an initial model of the disturbance. The disturbance model is 
then used to calculate the optimal state estimator. If the ini- 
tial LQG controller meets minimum stability requirements, 
then the controller is applied. The update procedure is im- 
plemented by updating the disturbance model to account for 
changes introduced by neighboring control loops. In other 
words, the disturbance model is updated to make the system 
more robust in frequency bands where modeling error is ex- 
cessive. Next, the LQG controller is redesigned using the 
new disturbance model and a more aggressive effort weight- 
ing model. The new controller is evaluated based on control 
effort and performance criteria. If the design is acceptable, 
then the online controller is updated. The process is then re- 
peated until the control effort approaches predefined limits, 
or until subsequent designs fail to improve control perfor- 
mance. Additional details pertaining to this approach can be 
found in Schiller [10]. 

6 Experimental Setup 

The control strategy was evaluated using real-time experi- 
ments performed on a flat aluminum skin-stringer panel. The 
1.17 m x 1.17 m x 0.00127 m panel, shown in Lig. 3, was 
segmented into 10 approximately equal size bays with nom- 
inal dimensions of 0.45 m x 0.19 m. Lor logistical reasons, 
transducers were only attached to the two bays shown in 
Lig. 4. The remaining bays were covered with lead- vinyl. 
Lor these tests, the stiffened panel was mounted in the trans- 
mission loss window in the Structural Acoustic Loads and 
Transmission (SALT) facility at NASA Langley Research 
Center. The panel was driven by a broadband acoustic dis- 
turbance, while radiated sound power was estimated using a 
discrete array of structural measurements. In particular, a 5 







Figure 3: Photograph of the stiffened aluminum test panel. 
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Figure 4: Diagram of the front of the panel. 


x 6 grid of miniature accelerometers was used to measure 
the structural response of the two bays. Acoustic power was 
then estimated using radiation modal expansion [11]. Data 
acquisition and control were implemented using xPC Target. 
Since the response was sampled at 3 kHz, 1 kHz anti-aliasing 
prefilters and smoothing postfilters were used on all inputs 
and outputs. 

The control transducers consisted of surface-mounted piezo- 
ceramic patches and accelerometers, as shown in Fig. 4. 
Specifically, 0.07 m by 0.29 m piezoelectric actuators were 
mounted in the center of each bay, and four accelerome- 
ters were located in a diamond pattern around each actuator. 
When integrated, the summed response from the accelerom- 
eters provided an estimate of the volume velocity of each 
bay. This transducer configuration was selected based on the 
controller/transducer complexity work performed by Gibbs 
etal. [12]. 


7 Results and Discussion 

Initially input/output data was acquired to characterize each 
subsystem. This was accomplished by driving the PZT ac- 
tuators with broadband random noise while recording the re- 
sponse of each sensor. The observer/Kalman filter identifica- 


tion (OKID) algorithm [13] was then used to compute 75 th 
order state-space models of each subsystem. Conservative 
initial LQG controllers were then designed and implemented 
on each subsystem. The controllers were updated indepen- 
dently using the procedure described in the previous section. 
In particular, 60 th order disturbance models were updated 
based on the closed-loop response, while the effort penalty 
was slowly reduced. This procedure was implemented itera- 
tively until the control system converged. 

The power spectra of the volume velocity estimates acquired 
on the upper bay are shown in Fig. 5. The open-loop response 
is shown with the thin blue line, the initial closed-loop per- 
formance is indicated with the thick red line, and the final 
closed-loop response is shown by the dotted black line. The 
poor initial performance at 145 Hz was due to the interac- 
tion between the local controllers. As the controllers were 
updated, the amplitude of the disturbance model increased 
around 145 Hz. This significantly improved the performance 
of subsequent designs. The final design reduced the 145 Hz 
peak by more than 10 dB and achieved a 4 dB integrated re- 
duction from 50 to 600 Hz. Since the same trends were also 
observed on the lower bay, those results are not shown here. 

The total radiated sound power from both bays is shown in 
Fig. 6. The final design achieved a 2.7 dB integrated re- 
duction from 50 to 600 Hz, as shown by the thick red line. 
Spillover above 700 Hz was due to the fact that the integrated 
and summed accelerometers provided a poor estimate of vol- 
ume velocity at high frequencies. When a more accurate es- 
timate of radiated sound power was used for control, the sys- 
tem achieved 12 dB peak reductions and a 3.6 dB integrated 
reduction in radiated sound power. In addition, the high fre- 
quency spillover from 600 to 800 Hz was reduced as shown 
by the dotted black line in Fig. 6. In this case, the control sys- 
tem had access to radiated sound power estimates generated 
in real time using the radiation modal expansion technique 
[ 11 ]. 

8 Concluding Remarks 

Decentralized controllers introduce unavoidable errors due 
to the unmodeled coupling between subsystems. In partic- 
ular, the modeling error can be expressed as the product of 
the complementary sensitivity function for the neighboring 
control loop multiplied by a term that quantifies the diago- 
nal dominance of the plant. While closed-loop performance 
can often be improved by penalizing control effort at the ze- 
ros of the local control models, that approach is not always 
sufficient. Therefore the concept of loop transfer recovery 
was used to develop a time-varying strategy. This iterative 
strategy was validated using real-time control experiments 
performed on a structural- acoustic system with poles close 
to the stability boundary, non-minimum phase zeros, and un- 
modeled dynamics. Experimental results showed that updat- 
ing the disturbance model based on closed-loop system mea- 
surements improved the robust stability and performance of 
the control system. In particular, the iterative control system 
achieved 12 dB peak reductions and a 3.6 dB integrated re- 
duction in radiated sound power from two adjacent aircraft- 
style bays. 





Figure 5: The power spectra of the measured response on the upper bay before control (thin blue line), using the initial 
controller (thick red line), and using the final design (dotted black line). 



Figure 6: The total radiated sound power from the two bays before control (thin blue line), using the standard control system 
with summed accelerometers (thick red line), and using a more accurate radiated sound power estimate (dotted black line). 


References 

[1] R.H. Cabell, D.E. Cox, G.P. Gibbs, ’’Interac- 
tion metrics for feedback control of sound radi- 
tion from stiffened panels”, Proceedings of the 44th 
AIAA/ASM E/AS C E/AH S Structures, Structural Dynam- 
ics, and Material Conference 7-10 April, Norfolk, Vir- 
ginia (2003) 

[2] G.P. Gibbs, R.H. Cabell, ’’Simultanteous active con- 
trol of turbulent boundary layer induced sound radia- 
tion from multiple aircraft panels”, Proceedings of the 
8th AIAA/CEAS Aero acoustics Conference 17-19 June, 
Breckenridge, Colorado (2002) 

[3] J.C. Doyle, ’’Guaranteed margins for LQG regulators”, 
IEEE Trans. Automatic Control AC-23: 756-757 (1978) 

[4] B.D.O. Anderson, J.B. Moore, Optimal Control: Lin- 
ear Quadratic Methods , (Prentice Hall, New Jersey, 
1990) 

[5] R.H. Cabell, M.A. Kegerise, D.E. Cox, G.P. Gibbs, 
’’Experimental feedback control of flow-induced cavity 
tones”, AIAA J. 44:1807-1815 (2006) 

[6] J.C. Doyle, G. Stein, ’’Robustness with observers”, 
IEEE Trans. Automatic Control AC-24:607-611 (1979) 

[7] R.R. Bitmead, M. Gevers, V. Wertz, Adaptive Opti- 
mal Control: The Thinking Man ’s GPC , (Prentice Hall, 
New Jersey, 1990) 

[8] G. Stein, M. Athans, ’’The LQG/LTR procedure for 
multivariable feedback control design”, IEEE Trans. 
Automatic Control AC-32: 105-1 14 (1987) 


[9] M. Morari, E. Zafiriou, Robust Process Control , (Pren- 
tice Hall, New Jersey, 1989) 

[10] N.H. Schiller, ’’Decentralized control of sound radia- 
tion from periodically stiffened panels”, PhD thesis, 
Virginia Tech, Blacksburg, Virginia (2007) 

[11] G.P. Gibbs, R.L. Clark, D.E. Cox, J.S. Vipperman, ’’Ra- 
diation modal expansion: Application to active struc- 
tural acoustic control”, J. Acoust. Soc. Am. 107:332-339 
( 2000 ) 

[12] G.P. Gibbs, R.H. Cabell, J. Juang, ’’Controller complex- 
ity for active control of turbulent boundary layer noise 
from panels”, AIAA J 42,1314-1320 (2004) 

[13] J. Juang, Applied System Identification , (Prentice Hall, 
New Jersey, 1994) 


